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ABSTRACT 



The lattice formulation of quantum gravity provides a natural framework in which non-perturbative prop- 
• i— i . 

erties of the ground state can be studied in detail. In this paper we investigate how the lattice results relate to 
the continuum semiclassical expansion about smooth manifolds. As an example we give an explicit form for 
the lattice ground state wave functional for semiclassical geometries. We then do a detailed comparison be- 
tween the more recent predictions from the lattice regularized theory, and results obtained in the continuum 
for the non-trivial ultraviolet fixed point of quantum gravity found using weak field and non-perturbative 
methods. In particular we focus on the derivative of the beta function at the fixed point and the related 
universal critical exponent v for gravitation. Based on recently available lattice and continuum results we 
assess the evidence for the presence of a massless spin two particle in the continuum limit of the strongly 
coupled lattice theory. Finally we compare the lattice prediction for the vacuum-polarization induced weak 
scale dependence of the gravitational coupling with recent calculations in the continuum, finding similar 
effects. 
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1 Introduction 



It is widely believed that an understanding of the properties of quantum gravitation would have important 
consequences in many areas of cosmology and high energy physics. Unfortunately approaches to quantum 
gravity based on linearized perturbation methods have had moderate success so far, as the underlying theory 
is known not to be perturbatively renormalizablc [1, 2]. A lack of perturbative renormalizability implies 
that an increasing number of counterterms needs to be added in order to make the theory finite order 
by order in perturbation theory. It has been recognized for some time though that the lack of perturbative 
renormalizability is not necessarily an obstacle in defining a consistent quantum theory [3] , as several simpler 
field theory models suggest [4] (most notably the non-linear sigma model above two dimensions) and recent 
rigorous results seem to support [5]. In the continuum non-perturbative renormalizability requires the 
existence of a non-trivial ultraviolet fixed point of the renormalization group. In the presence of a lattice 
momentum cutoff, the corresponding requirement is the existence of a phase transition with a divergent 
correlation length. 

In the simplicial lattice formulation of quantum gravity one proceeds in a way similar to ordinary lattice 
gauge theories, and introduces a lattice ultraviolet regulator which in principle allows for controlled, sys- 
tematic analytical [6, 7] and numerical [8, 9, 10] non-perturbative calculations of ground state properties, 
anomalous scaling dimensions and invariant correlations. Once the specific lattice action has been chosen, 
numerically exact results can be obtained on finite volume lattices which then need to be judiciously extrap- 
olated to the infinite volume limit using the well-established methods of finite size scaling. In this paper we 
address the basic issue of the relationship between recent lattice results and a variety of approximate pertur- 
bative and non-perturbative results obtained in the continuum formulation, with the intent of establishing 
a set of connections between the two formulations that go beyond weak coupling perturbation theory and 
the perturbative expansion about smooth manifolds. 

The starting point for a non-perturbative study of quantum gravity is usually a suitable definition of the 
discrete Feynman path integral. In the simplicial lattice approach one starts from the discretized Euclidean 
path integral for pure gravity, with the squared edge lengths taken as fundamental dynamical variables, 

^ = r Ii II dl l ®VU ex P I - E -kShA h +aS 2 h Al/V h + ...)[ (1.1) 

J ° s ij { h J 

(see reference [10] for notation). The above expression is supposed to represent a lattice discretization of the 
continuum Euclidean path integral for pure quantum gravity 

Zc = /n {Vai^y n cx p {-J^ x ^ ( a - \ r + 1 R ^ Rixvpa +■■■)} > (l2) 
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with k^ 1 = 8irG, G Newton's constant, and the • • • represent higher order curvature invariant terms. The 
Rcgge lattice action only propagates spin two degrees of freedom in the weak field limit, while the cosmological 
and measure terms contain only local volume contributions. The discrete gravitational measure in Zl can 
be considered as the lattice analog of the DeWitt [11, 12] continuum functional measure [13, 14, 15]. A 
cosmological constant term with bare A > is needed for convergence of the path integral [8, 16, 17]. 
Without loss of generality one can rescale the metric and set A = 1. The curvature squared terms (a — > 0) 
allow one to control the short distance fluctuations in the curvature, and as far as the functional measure 
parameter a is concerned most of the recent work has focused on the case a = (for more details on the 
choice of action and measure the reader is referred to the review [9] , and further references therein) . 

The present numerical evidence from the discrete model of Eq. (1.1) suggests that quantum gravity in 
four dimensions exhibits a phase transition in the coupling G between two physically distinct phases [10]: a 
strong coupling phase, in which the geometry becomes smooth at large scales, 

(9^) « (G>G C ) (1.3) 

with a vanishingly small average curvature in the vicinity of the critical point at G c , and a weak coupling 
phase 

(S M „> = (G < G c ) (1.4) 

in which the geometry becomes degenerate, bearing some resemblance to a dilute branched polymer. It is 
clear that based on its geometric properties, only the smooth phase is physically acceptable. The existence 
of a phase transition at finite coupling G is usually associated with the appearance of an ultraviolet fixed 
point of the renormalization group, and implies in principle non-trivial scaling properties for the coupling 
constant and invariant correlations. 

In the lattice theory the presence of a fixed point or phase transition is usually inferred (as in other lattice 
field theories) from the appearance of non-analytic terms in invariant local averages, such as the average 
scalar curvature 

<T f d 4 x . fri R(rr\ -> , , 

A n (k c - k) 4 "- 1 . (1.5) 



Without such singularities the lattice continuum limit cannot be taken, as one needs a divergent correlation 
length to define the lattice continuum limit. Indeed k c here is defined as the location of the non-analyticity 
in the partition function and its averages, the latter often obtained by differentiation with respect to a source 
or some other parameter. A precise determination of v then allows one to connect singularities in averages 
such as the one above to other long-distance properties of the theory. In particular the relation between 
the critical exponent v and the derivative of the renormalization group beta function at the fixed point 
fi'{G c ) = —1/u implies a scale dependence of Newton's constant (due to gravitational vacuum polarization 



effects) of the form 

G(r) = G(0) [l + c{r/O l/v + O((r/0 V n] , (1-6) 

where £ is a renormalization group invariant scale parameter and c a calculable numerical constant of order 
one. Detailed numerical studies of the Regge lattice gravity model give a value very close to v~ x = 3 [10]. 
Since one finds for the critical value G c w 0.626 in units of the lattice spacing, one would conclude that 
the lattice theory is not weakly coupled in the vicinity of the fixed point. It seems natural to interpret 
the momentum scale £ _1 as arising due to a gravitational analogue of dimensional transmutation, and it 
playing a role in gravitation similar to the universal scaling violation parameter A-p^ of QCD [18]. Other 
lattice approaches to quantum gravity based on discrete dynamical triangulations with fixed edge lengths 
and which give rise to a rather different phase structure are reviewed in [19]. 

In this paper we examine a number of fundamental issues that have a bearing on the relationship be- 
tween lattice and continuum models for quantum gravity. First we will consider the lattice analog of the 
semiclassical expansion for the ground state wave functional of continuum gravity. Within the continuum 
formulation, the semiclassical expansion about smooth manifolds with bounded quantum fluctuations allows 
one to exhibit in a clear and direct way the transverse traceless modes (or equivalently spin two modes) 
as the only physical gravitational degrees of freedom. Such an expansion is most easily carried out with 
the Euclidean functional integral approach, wherein the gravitational action is expanded in the weak field 
metric, and the resulting Gaussian integrals arc subsequently carried out. 

In trying to construct the lattice analog of the ground state functional for semiclassical gravity one has 
two options. The first procedure relies on constructing directly a lattice expression for the exponent of the 
ground state functional, obtained by transcribing the continuum expression in terms of lattice variables. 
The unique procedure we follow here is to proceed from the lattice expression for the gravitational action, 
specified on a fixed time slice, and supplemented by the appropriate vacuum gauge conditions. A crucial 
ingredient in this method is the correct identification of the correspondence between continuum degrees 
of freedom (the metric) and the lattice variables (the edge lengths). This correspondence is fixed by the 
relationship between the lattice Regge action and the continuum Einstein action, at least in the weak field 
limit. The resulting lattice expression is then equivalent to the continuum one by construction. 

The second procedure relies instead only on the expression for the lattice gravitational action, as computed 
in the weak field limit, and determines the explicit lattice form for the ground state functional for linearized 
gravity by performing explicitly the necessary lattice Gaussian functional integrals. The resulting discrete 
expression can then be compared to the continuum one by re-expressing the edge lengths in terms of the 
metric. It is encouraging that the resulting lattice expression completely agrees with what is found by using 
the previous method. 
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It is advantageous in performing the above calculation to introduce spin projection operators, which 
separate out the spin zero, spin one and spin two components of the gravitational action. As a by-product 
one can then show that the lattice gravitational action only propagates massless spin two (or transverse- 
traceless) degrees of freedom in the weak field limit, as is the case in the continuum. Furthermore, as 
expected, the lattice ground state functional for linearized gravity only contains those physical modes. 

In subsequent sections of the paper we examine systematically the relationship between recent non- 
perturbative results obtained in the lattice theory and corresponding calculations performed in the continuum 
theory. The latter suggest the presence of a non-trivial ultraviolet fixed point in G, and in some cases have 
even led to definite predictions for the universal critical exponent of quantum gravitation, which can therefore 
be compared quantitatively to the lattice results. 

Besides relying on the recent lattice and continuum results for quantum gravitation, one can also try 
to independently estimate the gravitational scaling dimensions using what is known based on exact and 
approximate renormalization group methods for spin zero (self-interacting scalar field in four dimensions) 
and spin one (Abelian non-compact gauge theories), for which a wealth of information is available on the 
critical indices. Based on this comparison, we will argue that these results too are consistent with what 
is known about the gravitational exponents in four dimensions. Finally we describe a simple geometric 
argument which interprets the value found for the gravitational exponent v~ x = 3. 



2 Ground State Wave Functional of Linearized Gravity 



According to the path integral prescription for Euclidean quantum gravity, the wave function of the state 
of minimum excitation for an asymptotically flat three-geometry with specified metric is 

* [ 3 5y -,t] =Af [ [dg^je- 1 ^ (2.1) 

where the integral is over all Euclidean four-geometries which are asymptotically flat and bounded at time 
t by an asymptotically flat hypersurface with induced metric 3 gij- 

Kuchaf [20] has given an expression for the ground-state wave functional of linearized gravity. In the 
vacuum or "Coulomb" gauge hik,k — 0, ha — 0, h i — the ground state functional is given by 

Mh% T ,t] = M cxp j-J^ J d 3 k Wk ft^(k) ^ T (k)J (2.2) 

where /i^ T (k) is a Fourier component of the transverse traceless part of the deviation of the three- metric 
from the flat three-metric in rectangular coordinates, 

Mx,t) = 3 9ij(x,t)-8ij , (2.3) 
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c<Jk = |k|, TV is a normalization factor, and lp = (167rG) 1//2 is the Planck length in a system of units where 
h = c = 1. Equivalently one can write, in real space and in terms of first derivatives of the fields, the 
expression 

M h7, t] = Wexp i^^~rjT J ^ J dV ■ (2-4) 

The above ground state wave function of linearized gravity was originally evaluated by Kuchaf using canonical 
methods [20] . Later the same formula was obtained by Hartle using the Euclidean path integral prescription 
[21] (sec Section 3). 

Electromagnetic Case 

It is instructive, in view of the calculations to be done for the gravitational case in the next sections, to 
consider as an aside the much simpler case of electromagnetism. Indeed a completely analogous set of results 
for the ground state functional holds for the electromagnetic case, and further brings out the relationship 
between the action and the quantity appearing in the exponent of the wave functional, which can be described 
as the time-slice action contribution. 

In the Coulomb gauge diAi = 0, A = the ground state functional is given by Kuchaf in terms of the 
transverse parts of the potentials only 

*o[A, t] = AAcxp J d 3 k k Af(k)AJ (k)) (2.5) 

or, equivalently, in terms of the B fields as 

*o[M = AAcxp (- JL / d 3 x J d s x , B Wj°jf ) . (2.6) 

It is easy to see that the expression in the exponent of Eq. (2.5) is related in a simple way to the original 
electromagnetic action. One has for the action appearing in the exponent of the Feynman path integral 

I[A„] = l -J d 4 xF^(x) F^(x) (2.7) 

which for a single mode reduces to 

h = \ (k^A v (k) - KA^k)) 2 . (2.8) 
In terms of the transverse projection of the field 

A J = (S l3 ) Aj (2.9) 
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one has for fco = Aq = 

k 2 Af (k)Af (k) = \ (kiAjQt) kjMk)) 2 (2.10) 
and therefore, after summing over all modes, 

J d^^Aj^AJ^) = J ±± (hA 3 (k)-k 3 Mk)) 2 (2.11) 

with Wk = Vk^- Thus the expression appearing in the exponent of Eq. (2.5) is the same as the time-slice 
contribution derived in Eq. (2.11). 

Spin Projections 



Returning to the case of the gravitational field one can follow the same procedure, first in the continuum 
and then on the lattice, and derive a lattice expression for the vacuum functional for linearized gravity. 
In the continuum the ground state functional for linearized gravity Eq. (2.2) can be obtained from the 
continuum action by suitably expanding the gravitational action in the weak field limit, and then imposing 
the appropriate gauge conditions. By later following the same procedure on the lattice, the corresponding 
discrete expression can be derived. 

The first step involves therefore the expansion of the continuum Lagrangian —^JgR in the weak field limit 

9nu = Vtiv + K hft U (2.12) 



with k = VSttG and \h^\ small. The quadratic part [2, 23] is given by 

£sym 2 dfih^^Oyhpp -\- — Op^h^vOphpv — dph^pdph^p — d^h^^Ophpp (2.13) 

with residual gauge invariance 

hpu -» V + d^ + d v ^ (2.14) 
and an arbitrary gauge function. For one mode with wave-vector k one has 

£sym ^ kphp V k v Hpp — kphp V kpYip V — kphupkpYivp — kph vlJ kphpp (2.15) 

= \ k l (-^23 - h 24 - h 34 + h 2 2h 3 3 + ^22^44 + h 33 h 4i ) -\ . (2.16) 

A gauge fixing term can be added of the form 

£fix = — 2 (k v hp_ v — —kphvis) (kph^p — — kphp P ) (2.17) 
giving for the combined gauge-fixed weak field action 

C-tot = £ S ym + £fix = — 9ph a/ 3dph a/ 3 + - d aa d ph p p . (2-18) 



However in the following we shall rely instead on the vacuum (or Coulomb) gauge fixing which gives the 
ground state functional for linearized gravity discussed previously. 

It will be advantageous to define three independent spin projection operators, which show explicitly the 
unique decomposition of the continuum gravitational action for linearized gravity into spin two (transverse- 
traceless) and spin zero (conformal mode) parts [24]. The spin-two projection operator is defined as 

2 1 1 

+ T^gk^kvkakfi + - {S^aSvp + 5^p5 va ) - -S^ v 5 a p , (2.19) 

the spin-one projection operator as 

P^ivctf) = {kfj.kaSv/3 + k^kpSva + k u k a S^p + k v kp8 ~^k^k v k a kp (2.20) 

and the spin-zero projection operator as 

1 11 

P ft°JaP = _ 3fc2 (k^KSap + kakpSpv) + -8^8 a p + T-^k^kvkakfi . (2-21) 

The sum of the three spin projection operators is then equal to unity 



As a result one can define for the metric three orthogonal fields of definite spin, the transversc-tracelcss (spin 
two) part 



1 



the longitudinal part (spin one) 



and the trace (spin zero) part 



such that their sum gives h 



h^ u — P^ahapPpv — -P^uPafjhfla (2.23) 



hT = V - P^ a h a php v (2.24) 



V = r^PcpKp (2.25) 



h = h TT + h L + h T . (2.26) 
Here we have defined the projection operator 

P^ = <J„„ - 1 d^d v (2.27) 



or equivalcntly in momentum space 

k^k u 

"fc2 



Ppu = <W - ( 2 - 28 ) 



8 



Using the three spin projection operators, the action for linearized gravity can then be re-written simply as 



r - -lh 



p(2) _ 2p (0) 



k 2 h a0 (2.29) 



Imposing the gauge condition hio = /i o = h>ik,k = and setting k — one obtains 

-\k* hjfh™ + \k* hj jh l (2.30) 

with the second (spin zero) vanishing after further imposing the trace condition ha = 0. The resulting 
expression is then identical, up to a factor, to the expression appearing in the exponent of the ground state 
functional of linearized gravity of Eq. (2.2). 

Lattice Transverse Traceless Modes 



In this section the lattice analog of the TT-mode action will be derived from the Regge lattice gravitational 
action, and a lattice expression for the gravitational wave functional will be given. As a first step one needs 
to perform the weak field expansion for the Regge action 

lR = ^2A h S h (2.31) 

h 

where Ah is the area of the hinge h, and Sh is the deficit angle at the same hinge. Following [7] each 
hypercube in a hypercubic lattice is divided up into 24 four-simplices, with vertices at (0, 0, 0, 0), (0, 0, 0, 1) 
... (1,1,1,1) (without loss of generality one can take the lattice spacing to be one). The lengths of the 
15 edges connecting the vertices i and j are denoted by kj, where i and j range from 1 to 15, with the 
coordinates of the endpoints interpreted as binary numbers (for more details, see Section 3). Next each link 
length is allowed to fluctuate by an amount 1 + e around the hypercubic lattice value. To lowest order in 
the edge fluctuation, the lattice action is given by a quadratic form 

J fl = ^X> Mtfei (2 ' 32) 

ij 

with M a local matrix connecting only nearest-neighbor points. In Fourier space one can write for each of 
the fifteen displacements e^' b '°' d , defined at the vertex of the cube with labels (a, b, c, d), 

e («,6,c,<0 = M>2) <> MCM <i e o (2 .33) 

with w\ = e lkl , uj 2 = e lk2 , ^4 = e lk:i , ^8 = e lk4 (we use the binary notation for u and e, but the regular 
notation for fcj. For one mode (one set of u's) one obtains therefore (see Appendix B in [7]) 

Qe\ + 16eg + I8ej + (W1W4 + ^2^4 + ^1^8 + ^2^8 + + W2W4 + Cuiujg + bj-iLO%)e\ei 
— (8 + 4w 2 + 4(I) 2 )eie3 — (2wi + 2u>i + 2a) 2 ^4 + 2w 2 tl>4)eie6 — (12 + 6^4 + 6a} 4 )e 3 e7 + • • • (2.34) 
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Each coefficient is real, as expected from the reality of the action. Thus, for example in the above expression, 
we have 



OJ1U4 + UJ2UJ4 + UJiUJs + W 2 0>8 + U)\U)4 + U! 2 L04 + Wl^8 + W2OJ8 

4cos(i(fc 3 -ki)) [cos(i(2fci + fc 3 + fc 4 )) + cos(±(2fc 2 + h + fc 4 ))] 
- 8 - 2kl - 2k\ - 2fcifc 3 - 2fc 2 fc 3 - 1k% - 2k x k 4 - 2k 2 k 4 - 2k\ + (9(fc 4 ) 



(2.35) 



To show the equivalence of the Regge action to the continuum Einstein action one needs to replace the e 
fields with metric components (or alternatively, as done in [7], use trace reversed metric components), with 
body principals expanded as 

ei = -1 + [1 + c^n] 1/2 
e 2 = -1 + [1 + 0J 2 h 22 ] 1/2 
e 4 = -1 + [1 + ^4/133] 1/2 
e 8 = -1 + [1 + LU 8 h M ] 1/2 , (2.36) 



face diagonals as 



e 3 = -1 + 
e 5 = -1 + 
eg = -1 + 
e 6 = -1 + 
-1 + 



1 + -WlW2(/lll + ^22) + /ll2) 



1 + -WlW4(/lll + ^33) + /ll3) 



1 + -UiUg(hu + ^44) + /114) 



1 + -W 2 W 4 (/l22 + /l33) + ^23) 



eio 



1 + -Lu 2 uj s (h 22 + ha) + /i 24 ) 



1/2 



1/2 



1/2 



1/2 



1/2 



ei2 



-1 



nl/2 



1 + -Ld 4 uj 8 (h 33 + ha) + h 34 ) 



(2.37) 



body diagonals as 
e 7 = -1 + 

en = -1 + 
eis = -1 + 
ei4 = -1 + 



1 + ^<x>ia>2W4(/lll + ft-22 + /l33) + ^((1 +^4)^12 + (1 + Wi)/l23 + (1 +^2)^13)) 
1 + ^1^2^8(^11 + ^-22 + h 44 ) + ^((1 + W 8 )ftl2 + (1 + Wi)/l 2 4 + (1 + U 2 )h 14: )) 
1 + ^WiW 2 W4(ftll + ft-33 + ^44) + ^((1 + Wg)/li 3 + (1 + Wi)/l 34 + (1 + W 4 )/li 4 )) 
1 + ^LU 2 LJ4LJs(h 22 + h 33 + /144) + ^((1 + UJ 8 )h 23 + (1 + W 2 )/l 3 4 + (1 + W4)/l24)) 



-,1/2 



-,1/2 



1/2 



1/2 



(2.38) 
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and finally hyperbody diagonals as 

r i 3 i 1/2 

ei 5 = -1+ 1 + -uiiu)2U4V8(hu + h 22 + h 33 + /144) + 7(^12 + /ll3 + hi4 + h 2 3 + h 2 4 + flu) ,(2.39) 



although the latter quantity is not needed, as it does not appear in the Regge action to lowest order in the 
weak field expansion. Each expression is then expanded out for weak h, giving for example 

ei = ^ wiftn + 0(h 2 ) 

e 3 = 7} h 12 + ^LJiuJ2(hii + h 2 2)+ 0(h 2 ) 

&1 = \ (hi2 + h 13 + h 23 ) + \ (wi/123 + W2/113 + W4/112) + 7, Luiuj 2 ^4(hu + h 22 + h 33 ) + 0(h 2 ) 
00 

(2.40) 

and so on for the other edges by permuting indices. Setting then uii — e lkl ... oj$ — e tki (we switch here from 
the binary notation for the w's to a normal notation for the fc's), the resulting answer is finally expanded 
out in k to give exactly the weak field expansion of the continuum Einstein action as given in Eq. (2.13) and 
Eq. (2.15), and completely parallels the procedure for recovering the continuum limit of the lattice action as 
described in [7, 26]. 

To derive a lattice expression for the ground state functional of linearized gravity one needs to compute 
the lattice gravitational action on a given time slice, and subsequently impose the appropriate discrete 
vacuum gauge conditions. This will then give the action contribution appearing in the exponent of the 
ground state functional for linearized gravity as it appears in Eq. (2.2). 

The first step involves therefore the imposition of the vacuum gauge conditions hik,k = 0, h on = h oi = 
which gives 

e 8 = 
1 

eg = 2 Ws ei 
1 

eio = 2 Ws e2 
1 

ei2 = 2 Ws e4 

en = \ (1 +uj s )e 3 - \ (1 - uj s )(iv 2 ei +wie 2 ) 
3 

e 13 = - (1 + w 8 )e 5 - - (1 - w 8 )(w 4 ei + wie 4 ) 

e 14 = \ (1 +uj s )e 6 - \ (1 - w 8 )(w 2 e 4 + w 4 e 2 ) (2.41) 
3 

and results in an action contribution of the form 

2e\ + 8el + (W1W4 + w 2 w 4 + ^1^4 + ^2 + ^4,)eie 2 - 2(lu 2 + lu 2 + 2)eie 3 

-2(wiW2 + &i&2 + W4 + <I> 4 )e 3 e4 + 4(w 2 + u> 2 + ^4 + ^4)e 3 e 5 H (2.42) 

11 



where the dots indicate again additional terms obtainable by permutation of indices. 

To verify that this is indeed the correct expression one can use the expansion of the e^'s in terms of the 
hij's, as given in Eq. (2.39), and then expand out the w's in powers of k. One obtains 

7}( k i h 23 - k\h 22 h 33 - 2k 1 k 2 h Y3 h 23 + 2k x k 2 h 12 h 33 + k\h\ 3 - k\h Y \h 33 + 2k x k 3 h 13 h 22 
—2k 1 k 3 hi 2 h 23 — 2k 2 k 3 h 12 h 13 + 2k 2 k 3 hnh 23 + k 3 h\ 2 — k 3 huh 22 ) (2.43) 

which can in turn be re-written as the sum of two parts, the first part being the transverse-traceless contri- 
bution 

\k 2 ^h(P 3 hP- l -P^{P 3 h))] = Jk 2 /^(k) /^(k) (2.44) 
W h V = Tr[ 3 /i(P 3 hP- ^PTr(P 3 h))} (2.45) 
with Pij = Sij — kikj/k 2 , and the second part arising due to the trace component of the metric 

-ik 2 Tr[PTr(P 3 /i)PTr(P 3 /i)] = k 2 /^-(k) fcjT.(k) (2.46) 

with h T = |PTr(P 3 h). In the vacuum gauge hik,k = 0, ha — 0, h M — one needs to further solve for the 
metric components h\ 2 , h\ 3l h 23 and h 33 in terms of the two remaining degrees of freedom, hu and h 22 , 

h 12 = - 1 - (h n k 2 + h 22 k\ + fenfc 2 + h 22 k\) 
2k\k 2 

h 13 = ~ 2fc 1 fc (hnk 2 - h 22 k\ - h\\k\ - h 22 kl) 
h 23 = - 1 (-hnk 2 + h 22 k\ - h xl k\ - h 22 k\) 

h 33 = -h n - h 22 ( 2.47) 

and show that the second (trace) part vanishes. For example, in terms of the e variables the vacuum gauge 
condition hik,k = then reads 

[2(1 - u>i) +w 2 (l - uj 2 ) +w 4 (l - UA)]ei+uJi{l-uj 2 )e2+ui{l-uJi)e4,-2(l-uj 2 )e 3 -2(l-w 4 )e< i = (2.48) 

and permutations. The above manipulations then show that the expression given in Eq. (2.42) is indeed 
the sought-after lattice analog for the continuum expression k 2 /i^ T hf^ appearing in the exponent of the 
ground state functional of lineraized gravity. 

We conclude this section by outlining an example of a potentially useful application for the above results. 
The explicit construction of the ground state wave functional of linearized lattice gravity in terms of lattice 
transverse-traceless modes makes it possible at least in principle to compare the lattice and continuum results 
in the limit of small curvatures, such as would be obtained for example from lattice simulations by imposing 
flat boundary conditions at spatial infinity. After imposing the boundary conditions by suitably restricting 
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the values for the edge lengths on the lattice boundary such that the deficit angle is zero there, one would then 
have to further enforce the lattice vacuum gauge conditions of Eq. (2.48) so as to finally make contact with 
the semiclassical lattice functional of Eq. (2.42). But no gauge fixing is required for determining invariant 
averages obtained via the partition function of Eq. (1.1), so in practice the gauge conditions would have to 
be imposed configuration by configuration, by progressively applying local gauge transformations [14] so as 
to gradually transform the edge lengths for each configuration to the lattice analog of the vacuum gauge. 
It is expected that after such a transformation the edge distributions on a fixed time slice should follow 
closely the distribution of Eq. (2.42), if indeed as expected the only surviving physical modes are transverse 
traceless. 

3 Ground State Wave Functional for Linearized Regge Calculus 

In the previous section the vacuum functional for linearized gravity was derived by evaluating the discrete 
Euclidean action on a fixed time slice, which was later supplemented by the appropriate gauge conditions. In 
this section the ground state functional will instead be derived by performing directly the discrete functional 
integration over the interior metric perturbations for a lattice with boundary. 

One of the main results of this paper is the analog of Hartle's continuum calculation of the ground state 
wave functional of linearized gravity, using Regge calculus. We shall now briefly outline Hartle's calculation 
since our calculation later is (intended to be) a discrete version of his. 

In linearized gravity, the Einstein action is expanded to quadratic order in deviations of the metric from 
its flat-space value. On a surface which becomes the flat surface t = constant when the metric perturbations 
are zero, we write the three-metric as 



and hij can be decomposed into a transverse traceless part, a longitudinal part and the trace. Since the 
physical degrees of freedom are the two independent components of hf^ , the transverse traceless part, the 
wave function on the surface can be written as 



(3.1) 



*o = *o[^(x), t]. 



(3.2) 



The Euclidean Einstein action is given by 
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l 2 P I(g) = - f d A x^g R - 2 f d 3 x^ K, (3.3) 

JM JdM 

and for linearised gravity, the Euclidean four-metric in the functional integral is written as 

g a p{x) = 6 a p + h af} {x) (3.4) 

and the action is expanded to quadratic order in h a p. The boundary dM is taken to be a flat slice in flat 
Euclidean space, and M is the region of flat Euclidean space to the past of this. The h a p are required to 
vanish at infinity so that g a p is asymptotically flat. 

The action is required to be gauge invariant, and gauge- fixing terms in the four-volume and on the surface 
are included in the functional integral. The metric perturbations are divided into conformal equivalence 
classes [22] by writing 

h a p{x) = <t> a p{x) + 2 5 a px{ri)- ( 3 - 5 ) 

The integration contour for \ is rotated to purely imaginary values to make the integral over \ converge. 
Then the field (j) a p is decomposed as 

4>af3 = <t> a p + fa/3, (3-6) 

where (j> a p is a solution of the linearised field equations which satisfies the gauge and boundary conditions. 
The unique solution is essentially that the spatial components are the hj? and the other components vanish. 
The integral over the f a p contributes only to the normalisation factor, as does that over \. The result for 
the ground state wave function of linearised gravity is 

MhJ 3 T , A = AAcxp (- J rf 3 k co k ^ T (k) ^ T (k)) , (3.7) 

where /i^ T (k) is the Fourier transform of /i^ T (x), ujk — |k|, M is a normalisation factor and lp — ^/T&ttG is 
the Planck length. This is exactly the same formula as that obtained by Kuchaf using canonical methods. 

Linearized Regge calculus can be implemented as the theory of the small fluctuations of edge lengths 
away from their flat-space values, in a tessellation of four-dimensional space using rectangular hypercubes 
subdivided into four-simplices. The methods of subdivision and the notation are described in detail in earlier 
work on linearised Regge calculus [7], the difference in this case being that we have a four-dimensional 
Euclidean space with a flat boundary. The binary notation in [7], which we shall also use here, comes from 
interpreting lattice vectors (0, 0, 0, 1), (0, 0, 1, 0), (0, 1, 0, 0) and (1, 0, 0, 0) as binary numbers, giving the Xi-, 
X2-, Xi- and xs-directions. For ease of notation here, we shall take the boundary surface to be xs — 0, and 
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the Euclidean four-space to be x% > (to avoid lots of minus signs). Unlike the continuum case, we shall 
take periodic boundary conditions in the 1-, 2- and 4-directions, while the space will be asymptotically flat 
in the 8-direction. With unit lattice spacing, the flat-space edge lengths will be 1,^2,^3 and 2, and the 
perturbed edge lengths will be written as in Section 2, as 

4 = 14(1 + 4), (3.8) 

where L is the flat-space edge length, e is a small perturbation, and in each case, the upper index j denotes the 
lattice point at which the edge is based and the lower index i denotes the direction of that edge (all in binary 
notation). (Note a small change in notation from [7] where the small perturbations were called 5.) Thus, for 
instance, the e's based at the origin and lying in the boundary hypersurface will be e°, e\, e§, e% e°, eg, e°. 

In brief, our method is to write down the action for the semi-infinite four-dimensional space and to 
perform a functional integral over the internal perturbations, leaving an expression in terms of the e's on the 
boundary. The aim is then to identify the quadratic expression in the surface e's with the discrete version of 
the integral of hf^ (k)h TTt: > (k). The calculation is long and tedious so we shall give relatively little detail, 
but enough for the reader to reproduce it if required! 

Interior Terms 

Consider an interior vertex, which for simplicity we shall label as if it were the origin. The e's for the edges 
based at this vertex will be e°, i = 1, 2, 15, and we first write the total quadratic contribution (the first 
non-vanishing order) to the action, which involves any of these e's. This will arise from the Regge action 
"SAhdfn for the hypercube based at the origin and from neighbouring hypcrcubes. This is given explicitly in 
Appendix B in [7]. 

The next step is to differentiate the action with respect to each of the e\ in turn to obtain their classical 
equations of motion. Below we give an example of an equation of motion of each type, for e\, e° and e° 
respectively: 

= 6e° - 4(e° + e° + eg) - 2(e 6 + e\ + e\ 2 ) - 4^ 2 + e 5 ~ 4 + e^ 8 ) 

+3(ef 2 + e n 2 + ef 4 + ef 3 4 + e n 8 + e^ 8 ) + e\ + e 3 8 + e\ + e 5 8 + e\ + e^ 

+e^ 6 + e7 6 - 2e^ + e^ 10 + e^ 10 - 2e w w + e7 12 + e^ 12 - 2e^ 2 (3.9) 

= 8e° - 2(e° + e°) - 3(e? + e^) + 2{-e\ + e\ + e\ ) + 2(-e\ + e 2 + e 2 ) 

-el - e 4 - e7 4 - + 2(e^ + + + e w *) 3(e7 4 + e n 8 ) (3.10) 
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= 6e° 7 - 2e% - 2e° 5 - 2e% + e\ + e\ - 2e\ + el + ej- 2e\ + e\ + e| - 2e\. (3.11) 

All other equations of motion may be obtained by cyclic permutations of the indices. 
We introduce new integration variables f- by 

<4=%+f'i, (3-12) 

where the e\ satisfy the equations of motion above. By subtracting each of the e\ times the corresponding 
classical equation of motion from the contribution to the action based at the origin, the e\ are completely 
eliminated, leaving only Gaussian integrals over the ff, which contribute only to the normalization. The 
same feature appears in the continuum. (Note that as in [7], cross terms of the form e^e J k , where j is 
a neighbouring lattice point of the origin, are assigned half to each of the lattice points involved.) This 
elimination of the interior terms would seem to hold independently of whether we impose periodic boundary 
conditions or asymptotic flatness. 

Boundary Terms 

Having integrated over all contributions to the action from interior vertices, we are now left with the contri- 
butions assigned to vertices on the x$ = boundary. These will consist not only of terms involving e's based 
at vertices on the boundary hypersurface, but also of contributions from vertices one layer in but assigned 
(as explained above) partly to the boundary layer. 

Suppose now that the origin is back on the boundary hypersurface. The total contribution to the action 
involving e° has more than 200 terms so will not be reproduced here. Clearly there are no terms involving 
with j — 8,9, 14 since the boundary is at x 8 = 0. All e\ terms with i or j = 8,9, 14 are as for 
interior vertices, and most of the other terms are half their interior values. 

Recall that the e's in the boundary three-surface (e°, i = 1, 2, 7) are fixed but e° with i = 8, 9, 14 are 
to be varied and have exactly the same classical equations of motion as before. We again simplify the action 
by subtracting e- 1 times its equation of motion, for each of these, and then eliminate e^e^ and e° 4 using 
their equations of motion. At the same time, we eliminate e° using a constraint identical to its equation 
of motion. This looks somewhat suspect, but the motivation and justification are as follows. In three- 
dimensional linearised Regge calculus performed in a manner completely analogous to the four-dimensional 
case in [7], the e° mode is not dynamical and satisfies a constraint which turns out to be identical to its 
equation of motion; this reduces the number of variables to six, the correct number. We apply this result to 
our three-dimensional boundary hypersurface. 
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As a first step in linking our position-space representation of the action to the momentum representation 
in the Kuchaf-Hartlc formula [20, 21], we take the Fourier transform in the 1-, 2- and 4-directions, which 
are those in which there are periodic boundary conditions. (In the 8-direction, we have contributions from 
only one other layer, the first interior one). The details are explained in [7] but there the complex nature of 
the Fourier transforms was not taken into account. Here our convention is that e°e° transforms to w a e°e°, 
where a = 1,2,4 and ui a — e lka . For consistency with this, the linear expressions in the e's in the equations 
of motion transform slightly differently from in [7], with ef transforming to w a e°. To simplify the notation, 
we immediately drop the tildes from the Fourier transforms and the superscripts 0. 

We have not yet eliminated es,eg,em and e\ 2 ; this is not straightforward in the way it was for en,ei3 
and ei4 as their equations of motion are simultaneous equations for the four e's. For example, the Fourier 
transforms of those for e$ and eg, with a, defined to be 1 — u%, are 

e8(2(|ai| 2 + |ct2 1 2 + 1 0:4 1 2 ) — oio 2 — o 2 oi — «ia 4 — 0401 — 0204 — 0:40:2) 

+2e 9 (-|o 2 | 2 - |o 4 | 2 + oio 2 + 0104) + 2eio(-|ai| 2 - |o 4 | 2 + o 2 oi + o 2 a 4 ) 
+2e 12 (-|a 1 | 2 - |o 2 | 2 + 040! + 040:2) 
= ei(2ai + 2a2 + 204 — 02 — 04 — 201O2 — 20104 + 01O2 + 0104) 
+e 2 (2o! + 2o 2 + 2o 4 — ot\ — a 4 — 2aia 2 — 2o 2 o 4 + o 2 Oi + 0204) 
+e4(2ai + 2o2 + 204 — ai — 02 — 20104 — 20204 + 0401 + 0402) 
+2e 3 (2o!0 2 - ai - o 2 ) + 2e 5 (2o!0 4 - ai - o 4 ) + 2e 6 (2a 2 o 4 - a 2 - o 4 ) 
+ef(— 2ai — o 2 — 04 + 2o 2 + 204) + e\{— 2a 2 — 01 — o 4 + 2oi + 204) 
+64(^204 - ai - a 2 + 2ai + 2a 2 ) 
+26^(0! + a 2 ) + 2e|(ai + a 4 ) + 2e|(a 2 + a 4 ), (3.13) 

es(— |o 2 | 2 - |a 4 | 2 + a 2 o:i + o 4 ai) + 2e 9 (|a 2 | 2 + |o 4 | 2 ) — 2eioa 2 ai — 2ei 2 o 4 Q:i 

= ei(a 2 + 04) + e 2 (2ai - a 2 - o 2 ai) + e 4 (2oi - o 4 - OiOi) + 2e 3 a 2 + 2e 5 04 

+e\{a 2 + 04) + e^(o 2 - 2ai) + e^(o 4 - 2oi) - 2e^a 2 - 2e^a 4 , (3.14) 

and those for eio and e\ 2 by cyclic permutations of the indices. 

These equations are not all independent. Adding the left-hand sides gives zero, while adding the right- 
hand sides gives 
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= 2(ei(|a 2 | 2 + |o 4 | 2 - «i«2 - Oi\Oii) + e 2 (|«i| 2 + |ct4| 2 - «i«2 - 0:20:4) 

+e 4 (|oi| 2 + |o 2 | 2 — a x a 4 - a 2 a 4 ) + 2e 3 a 1 a 2 + 2e 5 a 1 a 4 + 2e e a 2 a 4 ). (3.15) 

This expression is precisely the three-dimensional scalar curvature at the origin [27] and is constrained to be 
zero. (The constraint on a three-dimensional hypersurface also includes (trK) 2 and trK 2 terms but these 
terms are of higher order in the e's.) We can also show that a\ times the eg equation plus a 2 times the 
eio equation plus a 4 times the e\ 2 equation gives zero on the left-hand side but on the right-hand side, it 
gives half the difference between the three-dimensional curvature scalars at the origin and at vertex 8, and 
so again is zero. ( 3 R is also zero at the interior vertex 8 as the constraint is propagated into the bulk.) 

Thus we have a situation where only two of the equations arc independent, and where the consistency is 
guaranteed by the constraint. To solve the equations symmetrically, we take 

A = e 8 , (3.16) 



fi = aieg + a 2 e w + a 4 e 12 (3-17) 
to be arbitrary parameters. We then obtain 

eg = ^{B + \X 24 + 2a li i) 1 (3.18) 
e w = ^(C + \X 14 + 2a 2f i), (3.19) 

ei2 = ^(D + XX 12 + 2a 4 fi), (3.20) 

where S = |oi| 2 + 1 0:2 1 2 + |o;4 1 2 and Xij = \ai\ 2 + \aj\ 2 — cti&k — ctj&k with k — 1,2,4 and k 7^ B,C 
and D are the expressions on the right-hand sides of the eg,eio and e i2 equations. We now substitute for 
eg, eio and e\ 2 in the boundary action. The total coefficients obtained for A and p, are both multiples of 3 R 
at the origin and so vanish. 

What remains is an extremely long expression. We write the coefficient of ej in this as Wi + Ui (for 
"wanted" and "unwanted"!). We will now give the expressions for i = 1,3, all others being obtainable by 
cyclic permutations of the indices. 

4T,Wi = [ei(|a 2 | 2 + |ct 4 | 2 - 0:20:4 - 0402) 
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+e2 ( — 1«4 1 2 + 201O4 + 2a 2 a 4 — ol\6l2 — 0104 — 0402) 
+e 4 (— \a 2 \ 2 + 2aia 2 + 2a 2 a 4 — U\U 2 — 0104 — a 2 a 4 ) 
+2e 3 (a 4 o:2 - |ct4| 2 ) + 2e 5 (a 2 04 - | «2 1 2 ) + 2e 6 (oia 2 + oia 4 - 2a 2 04)](2 + S) 
— 2ei(|o 2 | 2 + |ct4| 2 — 0:20:4 — 0402) 
— 2e 2 (— |ct4 1 2 + 20104 + 20204 — ol\5l 2 — 0104 — 040:2) 
— 2e4(— |a 2 | 2 + 2aia 2 + 2a 2 a 4 — aia 2 — oia 4 — 0204) 
—463(0402 — |a 4 | 2 ) — 4ef (a 2 a4 — |a 2 | 2 ) — 4e|(oia 2 + 0104 — 2a 2 04), (3.21) 

2EW3 = [61(0204 - |a 4 | 2 ) + e 2 (ai04 - |a4| 2 ) 

+64(0104 + a 2 a 4 — 2oia 2 ) + 2e 3 |a 4 | 2 - 2^50204 — 2e 6 aia 4 ](2 + S) 
— 2ei(a204 — I04I 2 ) — 2e%(cticti — |a 4 | 2 ) 
—264(0104 + 0204 — 2axa 2 ) — 4eg|o4| 2 + \^a 2 ot 4 + 4egOi04, (3.22) 

SC/i = e 2 (|ai| 2 — 02|oi| 2 + 0104 + 04O2 — 010402 + a 2 ) 

+e 4 (|oi| 2 — o 4 |ai| 2 + oio 2 + o 2 o 4 — oio 2 a 4 + a 2 ) 
+2630102 + 2650104 + 2e60204 
— ef (0102 + 01O4 — 01O2 — 0104) 
+e 2 {— |oi| 2 + |a 2 | 2 — aia 2 + oia 2 — 0104 + 02O4) 

+ e4(— |ai| 2 + |ct4 | 2 — 0104 + 0104 — Oi02 + O2O4) 

-2e|aia 2 - 2e|aia 4 - 2e^a 2 a 4 , (3.23) 
EC/ 3 = 2(eia 2 oi + e 2 oi02 - e 4 aio 2 (l - a 4 ) + eiOio 2 + 6*0102 + e\ot\u 2 ). (3.24) 

We see that 

Et/i = ( 3 i?(0)- 3 i?(8))/2 

— ei(|o2| 2 + |o4| 2 — ct\u 2 — u\u 2 — 0104) 
+e 2 {a\ — |o 4 | 2 — oio 2 + oio 4 + a 4 a 2 + 0204 — 010402) 
+e4(o 2 — |o 2 | 2 — 0104 + oio 2 + o 2 04 — o 2 04 — 010204) 
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+ (e^ + e 8 2 + e^)(|o 2 | 2 + K| 2 - a x a 2 - aia 4 ). (3.25) 
Repeatedly using 3 i?(0) = 3 R(8) — 0, we then have the following expression for the boundary action: 

T,(eiUi + e 2 U 2 + e 3 U 3 + e A U A + e 5 U 5 + e 6 U 6 ) = 

(ef + 4 + el - (ei + e 2 + e 4 )) 3 i?(0)/2 + U, (3.26) 
where the remainder U is given by 

U = [(ei[e~i(aice2 — oi\oi 2 + ol\Ol 4 — aia: 4 ) 

+e 2 (|o 2 | 2 + a 2 — Oia 2 — oeide 2 + a 4 a 2 + 6t 2 a A — 0:10:40:2) 
+e 4 (|a4| 2 + a 2 — 01O4 — aia 4 + 02O4 + 02O4 — 010204)] 
+e 2 [ei(|ai| 2 4- a 2 — 0201 — oia 2 + a 4 ai + ai«4 — 020404) 
+§2(0102 — oio 2 + o 2 04 — 0:204) 
+§4(|a4| 2 + a 2 — a 2 a A — a 2 ct A + 01O4 + 0404 — 01020:4] 
+e4[ei(|ai| 2 + a 2 — 040:1 — 0104 + 0201 + 0402 — 020401) 
+e 2 (|o 2 | 2 + 03 — 0402 - a 2 a 4 + Oi0 2 + 0402 — 010402) 

+64(0104 — 01O4 + 02O4 — 0204)] 
+2e 3 [ei(o20i + oio 2 ) + e2(aia 2 + oia 2 ) + 64010204)] 
+2e5[ei(o40i + aio 4 ) + §2010402 + (34(0104 + 0104)] 

+2e 6 [§1020404 + §2(0402 + o 2 a 4 ) + §4(0204 + 02« 4 )] (3.27) 

The next step is to expand in powers of k, using = e lki , on — 1 — (in contrast to the previous 
section, here we keep the binary notation for the fcj's). For the remainder U above, we obtain 

U = 2i{kie 1 + fc 2 § 2 + fc 4 e 4 ) 3 i?(0) + 0(fc 4 ). (3.28) 

For W\ and W 3l we have 

2(kj + fcf + kl)Wi = (ei - ef)(fe| + k\ - 2k 2 k 4 ) 

-(e 2 - ef)(fc| + 3fcifc 4 + 3fc 2 fc 4 + hk 2 ) - (e 4 - e|)(fc| + 3fcifc 2 + 3fc 2 fc 4 + hk 4 ) 
+2(§ 3 - e|)(fc 2 fc 4 - k\) + 2(§ 5 - e|)(fc 2 fc 4 - fc 2 ) 
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+2(e 6 - 4)(kik 2 + hk 4 + 2k 2 k 4 ) + Oik 3 ), (3.29) 

{k\ + k 2 2 + kj)W 3 = (ei - ef){k 2 k 4 - k 2 4 ) + (e 2 - elXftifc* - *2) 

+ (e 4 - ef)(fcifc 4 + fc 2 fc 4 + 2fcifc 2 ) + 2(e 3 - e\)kl 

-2(e 5 - e^)fc 2 fc 4 - 2(e 6 - ef)fcifc 4 + 0(/c 3 ). (3.30) 
Thus our final expression for the boundary action is 

2(fc? + fc 2 2 + fc 4 2 ) [ei[(gl " ^ + " 2k2ki) 
-(e 2 - e 2 )(fc 4 + 3fcifc 4 + 3fc 2 fc 4 + k\k 2 ) - (e 4 - e 4 )(fc 2 + 3fcifc 2 + 3/c 2 /c 4 + k\ki) 
+2(e 3 - ef)(fc 2 fc 4 - kj) + 2(e 5 - ei)(fc 2 /c 4 - 
+2(e 6 - et)(fcifo + hk A + 2k 2 ki) + 0(k 3 )} 
+e 2 [...] + e 4 [...] 
+2e 3 [(e! - e?)(/c 2 fc 4 - kj) + (e 2 - e«)(fc 1 fc 4 - fcf) 
+ (e 4 - e\){kxki + k 2 k 4 + 2k x k 2 ) + 2(e 3 - e%)kl 
-2(e 5 - e|)/c 2 fc 4 - 2(e 6 - e^fc^ + 0(^ 3 )] 

+2e 6 [...] + 2e 6 [...]] (3.31) 

The coefficients of e 2 and e 4 can be obtained from those of e 1; and those of e 5 and e 6 from those of e 3 , by 
cyclic permutation of indices. 

This is to be compared with the expression for hf J T h TT ^ stated earlier, Eq. (2.42) supplemented by the 
gauge conditions of Eq. (2.48). Note that that expression does not distinguish between and e^-ej, and 
once this is taken into account, the expressions are identical (apart from an overall minus sign which arises 
from the fact that the Regge calculus expressions are calculated from the changes in dihedral angles rather 
than the deficit angles). Thus we can write our expression for the action as 

I = J d 3 k hl T {h TTi \Q) - F T «(8)) (3.32) 



or 



I = J d 3 k h^i-^h™*) = J d 3 k cu k hf J T h TTtj (3.33) 
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where u>k — &8 = V 'k\ + k\ + k\. This is using the behaviour of h(k) in a space with Euclidean signature, 
where, if the expansion is periodic in the 1-, 2-, and 4-directions, we shall have 

h(k) = C cxp \i(k\X\ + k 2 x 2 + k^x^) — k^xs], (3.34) 

with fcg = k\ + k 2 + k\ to satisfy the wave equation Oh = 0. 

The ground state wave function obtained is thus identical to the continuum result of Kuchaf and Hartlc. 



4 The Evidence for Spin Two 



In the weak field limit the lattice theory described by the partition function Eq. (1.1) is known to 
be equivalent to the continuum theory of a massless spin two particle, as embodied in Einstein's General 
Relativity with a cosmological constant term. One would hope that the local gauge invariance (continuous 
lattice diffeomorphism invariance) of the discrete gravitational action under metric deformations - taken 
sufficiently small so as not to violate the triangle inequalities [14] - would be powerful enough to ensure 
that the lattice theory still describes a regularized model for quantum gravity, even away from smooth 
manifolds. In this section we shall examine the evidence in support of the argument that the lattice theory, 
treated non-perturbatively and in the vicinity of the critical point at G c where the lattice continuum limit is 
formally defined, still describes a massless spin two particle. A comparison will be made between the lattice 
results and those obtained recently in the continuum using a variety of perturbative (2 + e expansion) and 
non-perturbative (renormalization group combined with a derivative expansion) methods. A second line of 
approach will be to compare the lattice results for the critical exponents of gravitation with what is known 
cither exactly or approximately for other spin values (0,1) in four dimensions, and look for a discernible 
trend. 

The cornerstone for these kind of arguments is the basic idea of universality. It is known that the long 
distance behavior of quantum field theories is to a great extent determined by the scaling behavior of the 
coupling constant under a change in momentum scale [3]. It is also well known that asymptotically free 
theories such as QCD lead to vanishing gauge couplings at short distances, while the opposite is true for 
QED and self-interacting scalar field theories in four dimensions. More generally, the fixed points of the 
renormalization group need not be at zero coupling, but can be located at some finite coupling, leading to 
non-trivial fixed points not necessarily accessible by perturbation theory, or possibly even more complex 
fixed lines and limit cycles [3, 28]. 
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The existence of non-trivial ultraviolet fixed points, at which the theory becomes scale invariant, corre- 
sponds in statistical mechanics language to the existence of one or more critical points. There the partition 
function exhibits non-analyticities and singularities caused by infrared divergences associated with a diver- 
gent correlation length. It has been shown [9, 10] that lattice gravity exhibits precisely such a transition 
where, for example, the curvature fluctuation 

diverges at some k c . Such a divergence signals a singularity in the partition function itself, since averages such 
as the average curvature 1Z and the curvature fluctuation xn ar e related to derivatives of Zl (of Eq. (1.1)), 
with respect to the gravitational coupling k = l/(16nG). 

Simple scaling arguments allow one to determine the scaling behavior of correlation functions from the 
critical exponents which characterize the singular behavior of local averages in the vicinity of the critical 
point. The appearance of a singularity in the free energy F(k) is caused by the divergence of the correlation 
length £, which close to the critical point at k c is assumed to behave as 

f = 1/m fc ~ A i {k c -k)~ v (4.2) 

and defines the exponent v. Since for the singular part of the free energy one expects F S i ng {k) ~ simply 
on dimensional grounds, one then obtains by differentiation with respect to k for the curvature fluctuation 

X n(k) ^ A XK (fc c -AO-* 2 -* • (4-3) 

The last expression allows, at least in principle, a direct determination of the critical exponent v. Large 
scale direct numerical studies of the lattice theory [10] give the value v = 0.335(9) for G c = 0.626(11), which 
suggests v = 1/3 for pure gravitation. 

Apart from a detailed comparison between critical exponents (which will be done later in this paper), 
a number of direct and indirect arguments can be given in support of the fact that the non-perturbative 
lattice theory still describes a massless spin two particle in the vicinity of the critical point at k c . Firstly 
the gravitational lattice action only propagates spin two (transverse traceless) degrees of freedom, as shown 
explicitly in the weak field expansion of the previous sections (the lattice gravitational functional measure 
is completely local, and does not contain any propagation terms, to any order of the weak field expansion). 
This result is further supported by rigorous work describing the convergence of the lattice action towards the 
continuum one for smooth enough manifolds [29, 30]. Secondly, the static interaction of two heavy particles 
of mass m described by two world lines kept at a fixed distance d has been shown to scale consistently as 
the mass squared, as expected for gravitational type interactions [31]. In the following we will explore this 
delicate issue further, by pursuing the connection with available non-perturbative results in the continuum. 
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Ultraviolet Fixed Point 



One can contrast and compare the lattice results with what one obtains for quantum gravity in the 
continuum. Since gravity is not perturbatively renormalizable in four dimensions, one has to go to a lower 
dimension (d = 2) where the perturbative expansion becomes meaningful, and expand about that dimension. 
Similar expansion have been shown to be quantitatively very successful in scalar field theories [36, 37], but the 
series are shorter and a significantly larger extrapolation is required in the gravitational case (for a general 
review of the diagrammatic field theory methods as applied to statistical mechanics models see [38, 39, 40]). 

In the 2 + e perturbative expansion for gravity [32] (earlier references can be found in [33, 34, 35]) 
one analytically continues in the spacetime dimension by using dimensional regularization, and applies 
perturbation theory about d = 2, where the theory is formally power counting renormalizable and Newton's 
constant is dimensionless. An expansion in the number of dimensions of course goes back to Wilson's 
original work [3], and since then similar methods have been shown to be quite successful in determining 
among others the critical properties of the 0(n)-symmetric nondinear sigma model above two dimensions 
[41]. This model is not perturbatively renormalizable either, yet describes a completely well-defined and 
physically relevant statistical spin system, namely the universality class of the 3-d Heisenberg ferromagnct. 
The same dimensional expansion methods have been extended with some success to fermionic models as well 
[42]. 

In the gravitational case the dimensionful bare coupling is written as Go = A 2 ~ d G, where G is di- 
mensionless and A is an ultraviolet cutoff (not to be confused here with the scaled cosmological constant), 
corresponding on the lattice to a momentum cutoff of the order of the inverse average lattice spacing, 
A <~ 1/Z - The method has of course its share of problems, as the Einstein action is a topological invariant in 
two dimensions, which leads to kinematic singularities in the propagator. In addition, to recover the physical 
case d=4 requires a rather bold extrapolation from two dimensions. The series themselves are rather short 
and strong assumptions need to be made about the nature of possible singularities in the complex coupling 
constant plane (in particular the absence of singularites close to d=3). Still, one can view this approach as 
providing some sort of gauge-invariant resummation of a specific set of subdiagrams which may or may not 
be ultimately relevant in d=4. 

A double expansion in G and e = d — 2 then leads above two dimensions to a non- vanishing beta function 

8G 

ftG) = = (d - 2) G - ft G 2 - ft G 3 + • • • , (4.4) 

and consequently a nontrivial ultraviolet fixed point in G, since ft > for pure gravity. Integrating Eq. (4.4) 
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close to the fixed point, one obtains for G > G c a non-perturbative, dynamically generated mass scale 

It should be noted at this point that Eq. (4.5) is essentially the same as Eq. (4.2), with slightly different 
notations. It also brings out the central importance of the exponent v, and how it relates to the scale 
dependence of the coupling G. The derivative of the beta function at the fixed point defines the critical 
exponent v (which to lowest order is in fact independent of (3 ), 

0{G C ) = -l/u . (4.6) 

In the previous expression m is an arbitrary integration constant, with dimensions of a mass, and has to be 
associated with some physical scale to be determined (as in QCD) by physical considerations (we will argue 
that it is the analog of Aj^-jt for gravitation). It would appear natural here to identify it with the inverse 
of a gravitational correlation length (£ = m _1 ), perhaps a length scale associated with some average long 
distance curvature (more on this later). The above renormalization group result also illustrates in a direct 
way how the lattice continuum limit should be taken. It corresponds to taking the ultraviolet cutoff A — > oo, 
and therefore G — > G Cl with m held constant. For a fixed lattice cutoff, the continuum limit is approached 
by tuning G to G c . 

The value of the universal critical exponent v has important physical consequences, as it directly deter- 
mines the running of the effective coupling G(fi), where ii is an arbitrary momentum scale. The renormal- 
ization group tells us that in general the effective coupling will grow or decrease with length scale r = 
depending on whether G > G c or G < G c , respectively. The physical mass parameter m is itself by defi- 
nition scale independent, and therefore obeys a Callan-Symanzik renormalization group equation, which in 
the immediate vicinity of the fixed point takes on the simple form 

ix — m( M ) = n { A m ll | G(p) - G c \ v } = (4.7) 

with A m a numerical constant. As a consequence, for G > G c , corresponding to the smooth phase, one 
expects for the running, effective gravitational coupling [10, 15] 

G(r) = G(0) [ 1 + c(r/0 1/v + O((r/0 2/ ") ] , (4-8) 

with c a calculable numerical constant of order one. 3 The physical renormalization group invariant mass 
m = determines the magnitude of scaling corrections, and separates the short distance, ultraviolet 
regime from the large distance, infrared region. As already mentioned in the introduction, there are in fact 

3 At very short distances r ~ lp one finds finite perturbative corrections to the potential as well, which can be computed 
analytically using weak coupling diagrammatic techniques [43]. 
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indications that in the Euclidean lattice theory only the smooth phase with G > G c exists (since spacetimc 
becomes branched-polymer like and therefore degenerate for G < G c ), which would then imply that the 
gravitational coupling can only increase with distance (this point will be discussed further in Section 5). 

Recently the continuum 2 + e expansion for gravitation has been pushed to two loops, giving close to two 
dimensions [32] 

2 2f) 

P(G) = (d - 2) G - -(25 - n f ) G 2 - (25 -n f )G* +■■■ , (4.9) 

for nf massless real scalar fields minimally coupled to gravity. After solving the equation (3(G C ) = to 
establish the location of the fixed point, one obtains for pure gravity (nj = 0) 

G c = (3/50)(d - 2) - (9/250) (d - 2) 2 + • • • (4.10) 

and therefore close to two dimensions 

v- 1 = -P'(G C ) = (d-2) + \(d-2f + 0(d-2) 3 , (4.11) 

which gives to lowest order v^ 1 = 2 independently of d, and = 4.4 at the next order in d = 4. The 
uncertainty in these results can perhaps best be judged by comparing to similar calculations in the scalar 
case, for which much longer series exist, and for which rather sophisticated resummation methods based 
on Pade-Borel transforms, conformal mappings, and incorporating asymptotic large order estimates, are 
available [36, 37] (the methods of statistical field theory are discussed in detail in [38, 39]). Unfortunately in 
general the convergence properties of the 2 + e expansion for the non- linear sigma model are not encouraging, 
even when comparing to well-established results in d = 3 (e = 1) [41]. 

The 2 + e expansion is not the only method that has been applied in the continuum to extract quantitative 
informations about non-perturbative properties of gravitation. In this context we should mention another 
set of related results for the critical exponents of quantum gravitation. Recently in a separate, approximate 
renormalization group calculation based on the Einstein- Hilbcrt action truncation [46] one finds in the limit 
of vanishing bare cosmological constant v~ x = 2d(d — 2)/(d + 2) = 2.667 in d = 4, and v~ x w 1.667 in 
a more elaborate truncation. In this paper the sensitivity of the results to the choice of gauge fixing term 
and to the specific shape of the momentum cutoff is investigated as well. These more recent results extend 
earlier calculations for the exponent v done by similar operator truncation methods, and described in detail 
in references [44, 45]. A quantitative comparison of these various continuum results with the lattice answer 
for v~ x will be postponed until later in this paper. 

Geometric Argument for v — 1/3 
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A simple geometric argument can be given in support of the exact value v = 1/3 for pure quantum 
gravitation. The vacuum polarization induced scale dependence of the gravitational coupling G(r) as given 
in Eq. (4.8) implies the following quantum corrected static gravitational potential, for a point source of mass 
M located at the origin, 

V(r) = - G(r) — = - G(0) — \l + c (r/^ + 0({r/t) 2 ' v ) 1 (4.12) 
r r L J 

and for intermediate distances l p -C r <C £. As a result, the vacuum polarization effects due to virtual 

graviton loops cause an effective anti-screening of the primary gravitational source M. Thus the effect of the 

running gravitational coupling G(r) is to give rise to a new non-perturbative quantum contribution to the 

potential, proportional to r 1 /" -1 . Remarkably for v = 1/3 the additional contribution, now proportional to 

r 2 , can be interpreted as being due to what ultimately appears as a uniform mass distribution surrounding the 

original source. Its origin lies with a non-perturbative graviton vacuum polarization contribution, localized 

around the point source, and of strength 



P° = TZc3 ■ (4-13) 



3cM 
4<3 

Of course such a simple geometric interpretation fails unless the critical exponent v for gravitation is exactly 
one third. In fact in any dimensions d > 4 one would expect based on the geometric argument that 
—p'(G c ) = v~ x = d — 1, if the leading correction to the gravitational potential is due to a uniformly 
distributed, anti-screening cloud of virtual gravitons. These arguments rely of course on the lowest order 
result V(r) <~ J d d ~ 1 p e lp ' x /p 2 ~ r 3 ~ d for single graviton exchange in d > 3 dimensions. 

Equivalently, the running of G can be characterized as being due to a tiny non-vanishing (and positive) 
non-perturbative gravitational vacuum contribution to the cosmological constant, with 

k = -|s- ( 4 - 14 ) 

and therefore an associated effective curvature of magnitude 1Z ~ GA^ <~ GM/^ 3 . It is amusing that for a 
very large mass distribution, the above expression for the curvature can only be reconciled with the naive 
dimensional estimate 1Z ~ l/£ 2 , provided for the gravitational coupling G itself one has G ~ £/M [47]. 

Random Gravitational Paths 



Within the Feynman path integral formulation of quantum field theory, a well known relationship exists 
between the properties of random paths and those of field correlations (see for example [48]). In this section 
the analogy will be exploited in trying to gain more insight on the specific values for the gravitational critical 
exponents. 
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In the simpler case of self-interacting scalar field theories a rigorous argument can be given [50] based 
on an exact geometric characterization of criticality in the A</> theory and the Ising model, in and above 
four dimensions. The key element of the argument lies in the recognition of the fact that random walks 
representing the propagation of free particles in Euclidean space-time have fractal [49] dimensions dn = 
v^ 1 = 2, with vanishing probability of self-intersection above d = 4. As a consequence these models are 
governed by mean field theory above d = 4, with mild logarithmic corrections to free field behavior at d = 4. 
They provide rigorous support for the original claim that in the infinite cutoff limit all scalar field theories 
are trivial in four dimensions [3] . 

Let us first illustrate these results for the simplest case of a free scalar field in d dimensions with action 

S = \ f <j>(x)M(x,y)<t>(y) . (4.15) 

J x,y 

On a lattice one has for the matrix = — S^, where S is the (nearest-neighbor) hopping part, and the 
rest is the diagonal part Dij — (2d + m 2 ) Sij, with d the dimensions and mo the bare mass. The propagator 
connecting point 1 to point 2 is then given in terms of the kernel S by 

1 °° S 

Gl2 = ^TYd^^TTd (416) 

u n=0 u 

or, equivalently, in terms of a sum over paths 

da = J2 e- m/l2 ( path > ~ r^ 2 (d " 1)/2 e~ ri2 ^ (4.17) 

paths 1— »2 

where ^ 2 (path) is the length of the random path connecting points 1 and 2. In the second part of the 

expression we have indicated the asymptotic behavior of the free propagator for large distances, which 
brings in the correlation length £ = m^ 1 . 

In its simplest form, the lattice partition function needed to generate the above random curves is given 

by 

Z(J3) = AfJ Tl dDXi cxp|-/3^|X,-X 4+1 r| (4.18) 

which for a = 1 generates closed (Xjv+i = Xi) piecewise linear curves embedded in D Euclidean dimensions. 
For a = 2 it is equivalent to the generating function for a one-dimensional D-component massless field theory 
with unit lattice spacing, with infrared divergences appearing as the size N goes to infinity. In the limit of 
a large number of steps N one obtains for the size of the random walk 

1 N 

<X 2 >=— V<X^> - N 2/dH (4.19) 

i=l 

with dn = 2 for free Brownian motion, independently of a. The Hausdorff dimension du characterizes the 
deviation of an ensemble of random paths from what one would expect based on their topological dimension 
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of one. Furthermore dn is known to be a universal number, i.e. independent of the specific choice for the 
measure over random paths, and describing general geometric properties of random curves in the limit of 
very long paths. The relation v^ 1 = dn for free fields follows as a direct consequence of the representation 
of the field correlation function in terms of a sum over random paths with fixed endpoints, as given by 
Eq. (4.17). 

Below four dimensions non-trivial continuum behavior is expected for scalar fields, including non-Gaussian 
exponents and non-trivial fractal dimensions. Thus for example for an interacting scalar field in two dimen- 
sions v = 1 (Onsager solution of the Ising model), while for a self-avoiding random walk one has instead 
da = v~ x = 4/3 exactly in d = 2 [53]. Other non-trivial constraints, such as the requirement that the 
random walks do not back-track, are also expected to change the fractal dimension dn- 

In the gravitational case one is dealing with random paths associated with a massless particle of spin 
two. As a result new constraints on the nature of the random paths come into play, which are not present in 
the simpler case of a spinless scalar field. As discussed in Fcynman and Hibbs [54], these constraints are in 
fact already rather complicated for the case of a particle of spin one-half, and give rise even in this simplest 
case to a set of nontrivial complex weights needed to correctly reproduce the continuum expression for the 
Dirac propagator. In four dimensions such paths involve Dirac projection operators 1 ± 7^ [55]. 

On general grounds one would then be inclined to identify the value 1/ = 1/3 found for four-dimensional 
gravitation with a fractal dimension of random gravitational paths dn = 3. Unfortunately (or fortunately) 
the value u = 1/3 itself does not correspond to any known field theory or statistical mechanics model in four 
dimensions. For dilute branched polymers it is known that v — 1/2 in three dimensions [56], and v = 1/4 at 
the upper critical dimension d = 8 [57], so one would expect a value close to 1/3 somewhere in between. A 
value for the fractal dimension close to one would indicate the paths have almost linear Euclidean geometry, 
while at the opposite end a very large fractal dimension would indicate the paths are largely collapsed to 
a very small region about the origin. The paths in this latter case are highly folded and to some extent 
self-intersecting. Therefore the value dn = 3 found for quantum gravitation would suggest a far greater 
degree of folding compared to the spinless case, for which dn = 2. 

One could further develop these arguments and, in analogy with the scalar case, conclude that below 
six space-time dimensions two random gravitational paths will have a non-vanishing probability of self- 
intersection. These arguments would imply that the "upper critical dimension" for gravity is six, above 
which the theory becomes in some sense non-interacting and therefore trivial. 4 Unfortunately this argument 

4 One can obtain a direct estimate for the upper critical dimension by equating 2v~ 1 (d) = d, which after interpolating in d 
the known Rcgge lattice results gives as a solution the surprisingly low value d ~ 2.929 (see also Fig.l). This in turn would lead 
to the somewhat paradoxical conclusion that quantum gravity, in spite of being perturbatively non-renormalizablc, is weakly 
interacting in the infrared above three dimensions. In the same sense that self-interacting scalar field theories, in spite of not 
being perturbatively renormalizable, become weakly interacting at low energies above four dimensions. 
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is probably flawed, as dn = 3 holds only in d = 4, and presumably not at the upper critical dimension, if one 
indeed exists. On the other hand, the large-<i geometric estimate discussed previously, v^ 1 = d— 1, equates 
to two (the fractal dimension for an unconstrained, spinless random walk) in three space-time dimensions, 
where it is in fact known that there can be no propagating, genuine spin two degrees of freedom. 

Approaching Quantized Gravitation from Spin Zero and Spin One 

While only limited results exist for the non-perturbative scaling dimensions of quantum gravitation in 
four dimensions, the same is not quite true for spin one (compact Abelian gauge theory) and spin zero (self- 
interacting scalar field theory). It has been known since the work of Wilson [3] that all local one-component 
scalar field theories in four dimensions are described by the Gaussian fixed point. The fact that these theories 
become non-interacting (up to logarithmic corrections) at large distances in four dimensions implies that the 
critical exponents and scaling dimensions coincide with those of a free field. In particular one finds for the 
universal critical exponent 

v- 1 =2 (s = 0) (4.20) 

a result which in fact can be proven rigorously [50] (for some unresolved issues and an unconventional point 
of view regarding the self- interacting scalar field theory in four dimensions see [52]). In three dimensions 
the interacting A</> 4 scalar field theory shares the same universal long distance properties with the Ising 
ferromagnet. For both incarnations a wealth of numerical and analytical data exists on the critical exponents, 
and for the purposes of the present discussion the relevant result here is v~ x w 27/17 = 1.5882 [37]. 

In the massless spin-one case the results for the critical exponent v are somewhat less unambiguous, and 
furthermore no exact results are available yet. An analytical variational real space renormalization group 
analysis of the Abelian U(l) lattice gauge theory in d = 4 [59] gave 

v- 1 = 2.496(7) (s = 1) . (4.21) 

The errors there can be estimated from an analysis of the results for v using the same renormalization group 
methods in the 3-d U(l) spin system, which gave v = 0.6702(6) [59], compared to the present best theoretical 
value [37] v = 0.6698(15) based on the e-expansion about four dimensions as well as the 3-d </> 4 field theory, 
and also in good agreement with the latest experimental value v = 0.6706(5), as quoted again in the recent 
comprehensive review [37]. 

More recently [60] detailed numerical simulation studies of the Abelian U(l) compact lattice gauge theory 
have been performed with various additional action terms, besides the standard plaquette term. Away from 
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Reference 


v 1 in d = 3 


v~ x in d = 4 


HW93 [26] 


1.67(6) 


3.34 


H81 [59] 




2.991 


JLN96 [60] 




3.05(13) 


CF96 [61] 




3.48(12) 


exact 


1.5882 


3 



Table I: Critical Exponent u^ 1 for a massless spin two particle in four dimensions, as obtained indirectly either by 
extrapolation from other dimensions (d=3 in row 1) or from information on other spin values (rows 2-4)- Included 
in the table is also one direct (lattice) determination in d = 3. The un-weighted average of all extrapolated values 
listed in the second column is v~ x = 3.22. 

the "Wilson line" (where the adjoint coupling 7^0, and where the transition, being first order, is more 
difficult to analyze regarding the true singularity of the free energy located at the end of the metastable 
phase) the value v~ x = 2.53(6) is found for an adjoint gauge coupling 7 = —1/2. Closer to the first-order 
Wilson line 7 = they find the value h>~ x = 2.74(6), but which should really be discarded in view of the first 
order nature of the intervening transition for 7 = [62], unless as mentioned above a more refined analytic 
continuation towards the true critical point is performed, in order to extract the required critical exponent 
characteristic of the end-point singularity. 5 

Given the above values for the critical exponent v~ x for the massless spin zero and spin one fields, it is 
tempting to use them to estimate independently the scaling dimensions for gravitation, using 

v~\s) = svj^ + (1-aK-Jo . (4.22) 

(see Table I). With v~ x = 2 for s = and v~ x = 5/2 for s = 1 one then obtains v~ x — 3 for spin two, in 
good agreement with the previous discussion. In addition the simple formula 

v~ x = 2+| (4.23) 

gives for the exponent a/v = (2 — dv)jv = s in four dimensions, and therefore a divergence of the second 
derivative of the free energy of the remarkably simple form C ~ £ s . 

Yet another, independent way of estimating the critical exponent for four-dimensional quantum gravi- 
tation involves looking at the two lower dimensional cases of pure gravity in d = 2 (where v~ l = 0) and 
pure gravity in d = 3 (where v~ l w 1.67) [26]. A linear extrapolation to four dimensions would then gives 
v~ x = 3.3 which is quite consistent with what has been said in the previous discussion. It is worth noting 

5 One might wonder to what extent the quoted numerical results for the critical exponents, often obtained on relatively small 
lattices, are reliable. To further estimate the uncertainties in the numerical determination for v in d = 4 one can for example 
compare to a recent high-accuracy determination of v in the spin-zero case (Ising model) in d = 4, which yields v~ 1 = 1.992(6) 
[58] on comparable lattice sizes, and which should be compared to the expected exact value u^ 1 =2. 
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Reference 


v 1 in d = 3 


v 1 in d = 4 


HW93 [26] 


1.67(6) 


- 


H92 [10] 




3.08(62) 


H00 [10] 




2.98(7) 


AK96 [32] 


1.6 


4.4 


RS02 [44] 


1.11(5) 


1.68(26) 


RL02 [45] 




2.8(6) 


Li03 [46] 


1.2 


2.666 


exact 


1.5882 


3 



Table II: Direct determinations of the critical exponent v~ x for quantum gravitation, using a variety of analytical and 
numerical methods in three and four space-time dimensions. The un-weighted average of all direct determinations for 
quantum gravitation in four dimensions listed above gives v^ 1 — 2.93. 

here that the value for the exponent for three-dimensional gravity is tantalizingly close to the scalar field 
case. In the 2 + e-expansion one finds v~ x = 1.6 while some relatively old direct numerical simulations in 
d = 3 give v~ x = 1.67. Both values are quite close to the 3 — d scalar field exponent v^ 1 = 27/17 = 1.5882 
[37], which would be in line with the conjecture that in the infrared limit three-dimensional gravity belongs 
to the same universality class as a self-interacting single-component scalar field, with the scalar curvature 
playing the role of the scalar field R <~ U<p, as in fact suggested some time ago by the authors of reference 
[63]. 

As one last exercise one can look at the case of fractional spin, which presumably corresponds to massless 
self-interacting fermions. In the spin one half case, which should apply to self-interacting fermions in four 
dimensions (such as those represented by the non-renormalizable 4-d Gross-Neveu [42] and similar four- 
fermion models), one obtains v~ x = 9/4 = 2.25 in d = 4, which should be compared to the known values 
v^ 1 s=a 27/17 = 1.5882 in d = 3 (interacting 3-d fermions as described by Ising model exponents), and v = 1 
in d = 2 (based on the rigorous equivalence between the two-dimensional critical Ising model and a free 
Majorana fermion). Had one extrapolated linearly these known results to four dimensions, one would have 
estimated 2.18, a value quite close to 9/4 (to within three percent). It is of course not obvious at this point 
how to interpret the above result in terms of a fermion random walk, which would have fractal dimension 
dh = v~ x = 9/4. But the trend in the exponents is at least consistent with the expectation that the fractal 
dimension increases with embedding dimension d, as there are more dimensions to expand into. 

The various estimates for the critical exponent are compared in Tables I (indirect determinations) and 
II (direct determinations). Table II provides a list of critical exponents for gravitation as obtained by 
direct perturbative and non-perturbative methods in three and four dimensions. As mentioned before, 
direct numerical simulations for the lattice model of Eq. (1.1) in four dimensions give for the critical point 
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Gravitational Exponent v 1 
T 




Fig. 1 . Gravitational critical exponent v~ x = —fi'(G c ) as a function of dimension. Direct determinations from 
the Regge lattice (small circles at two, three and four dimensions), in the continuum using renormalization group 
truncation methods (squares), and by extrapolating lattice results from lower spin (triangles) are compared (see 
Tables I and II). The solid line is an interpolation through the Regge lattice results, incorporating the asymptotic 
behavior d — 1 for large d. The thin-dotted line is the analytic 2 + e result of Eq. (4.11). The dotted line is the 
continuum renormalization group result of [46] . The origin, methodology and comparison of the these various results 
is discussed further in the text. 



v- 1 


d= 2 


d = 3 


d = A 


spin s=0 


1 


1.588 


2 


spin s=l 








2.5 


spin s=2 





1.588 


3 



Table III: Summary table for the critical exponent v 1 as a function of spin and dimension. 
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G c = 0.626(11) in units of the ultraviolet cutoff, and one finds [10] 

v- 1 = 2.99(8) (s = 2) (4.24) 

which is used for comparison in Table II. The fact that the critical coupling G c is not small shows incidentally 
that the lattice theory is not weakly coupled close to the transition point. 6 

To conclude this section, one can reverse the line of the above arguments relating to the critical exponents 
for gravitation, and instead estimate the spin of the massless lattice graviton by considering the dependence 
of the measured exponent v on the spin. Assuming a linear dependence of the exponent v~ x on the spin and 
using the most accurate values at s = and s = 1 one obtains, from v^ 1 = 2.98(7) [10], about s = 1.98(3), 
which is quite close to the expected value of spin two. 



^Furthermore, the critical point obtained from the analytic continuation of the strong coupling (small k) branch of the free 
energy lies at the end of the metastablc phase of the Euclidean theory, which is not necessarily a concern here as one is ultimately 
interested in the pscudo-Ricmanniann theory. Indeed one would not have expected otherwise, in view of the well-known and 
seemingly un-avoidable conformal instability of the Euclidean theory. 
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5 Exponents and Long-Distance Quantized Gravitation 



The result of Eq. (4.8) implies that the gravitational constant is no longer constant as in the classical 
theory, but instead slowly changes with scale due to the presence of weak vacuum polarization effects, 

G(r) = G(0) [ 1 + c(r/0 1/v + O((r/0 V n ] , (5.1) 

The exponent v, related to the derivative of the beta function evaluated at the non-trivial ultraviolet fixed 
point via the relation /3'(G C ) = = —3 (see the previous discussion in Section 4 and reference [10]) 

for pure quantum gravitation, is supposed to characterize the universal long-distance properties of quantum 
gravitation, and is therefore expected to be independent of the specifics related to the nature of the ultraviolet 
regulator, introduced to make the quantum theory well defined. 

The mass scale m = in Eq. (5.1) determines the magnitude of quantum deviations from the classical 
theory, and separates the short distance, ultraviolet regime with characteristic momentum scale /i <C m 
where non-perturbative quantum corrections arc negligible, from the long distance regime where quantum 
corrections are significant. It should be emphasized here that most of these considerations are in fact quite 
general, to the extent that they rely mainly on rather general principles of the renormalization group and 
are in fact not tied to any particular value for the exponent v, although the value v = 1/3 clearly has some 
aesthetic appeal. Furthermore the dimensionless constant c is, at least in principle, a calculable number. In 
[64] c was estimated from the curvature correlation function at c = 0.014(4), while more recently in [31] it 
was estimated to be c — 0.056(27) from the correlation of Wilson lines. It is important to note that while 
the exponent v is universal, c in general depends on the specific choice of regularization scheme (i.e. lattice 
regularization versus dimensional regularization or momentum subtraction scheme). 

It is worthwhile to note that the result of Eq. (5.1), which as discussed in Section 4 is a direct consequence 
of Eq. (4.3) and the value for v found in the lattice theory (defined by the partition function of Eq. (1.1) with 
higher derivative coupling (a — > 0) and functional measure parameter a — 0), only applies to the simplest 
case of pure Einstein gravity with a bare cosmological term. 7 But one does not expect this to be the correct 
theory at sufficiently short distances r ~ lp, where higher derivative curvature terms will come into play, 
either through direct inclusion and gravitational radiative corrections, or via matter field and the conformal 
anomaly. In this limit the gravitational potential will be further modified by exponential and logarithmic 
terms, as discussed in reference [65]. 

Let us recall here that in SU(N) gauge theories and in particular in QCD, the theory of the strong 
interactions, a similar set of results is known to hold [18]. The crucial difference lies in the fact that there 

Tight matter fields will modify the exponent u, and therefore the result of Eq. (5.f), provided their mass is small enough 
to contribute significantly to vacuum polarization loops, m ~ 
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the scale evolution of the coupling constant can be systematically computed in perturbation theory due to 
asymptotic freedom, a statement which reflects the fact that such theories become free at short distances 
(up to logarithmic corrections). In non-Abelian gauge theories one has for weak coupling theory 

1 + 2/? log (JL.) + ... (5.2) 



9 2 W 9 2 (^ms) " b U M/ 
with fa = YjT^-(-4p — where TV is the number of colors, n/ is the number of flavors of massless fermions, 
/j, = 1/r is an arbitrary momentum scale, and A-p^ = 200MeV is a scale parameter which determines the 
size of scaling violations. The ■ ■ -'s indicate higher order loop corrections. Of course QCD does not determine 
Ajgj (it appears as an integration constant of the Callan-Symanzik renormalization group equations), and 
therefore it has to be fixed by experiment from a measurement of the size of scaling violations, i.e. via the 
observed deviations from free field behavior at sufficiently high energies. It is a remarkable fact that a good 
fraction of QCD and electroweak standard model phenomenology simply follows from the result in Eq. (5.2) 
and its electroweak analog. 8 

If one pursues in a straightforward way the analogy with non-Abelian gauge theories one is led to conclude 
that in quantum gravitation the quantity £ plays the same role as A-p^ in QCD, £ <-> Ajjg. One major 
difference between the two theories lies of course in the fact that in one case the ultraviolet fixed point is at 
the origin g 2 = 0, while in the other it is not. As a result one has logarithmic quantum corrections to free 
field behavior in QCD, but power law corrections in gravitation. 

To determine the actual physical value for the non-perturbative scale £ further physical input is needed. 
It seems natural to identify l/£ 2 with either some average spatial curvature, or perhaps more appropriately 
with the Hubble constant determining the macroscopic expansion rate of the present universe [15, 10], via 
the correspondence 

£ = , (5-3) 

in a system of units for which the speed of light is equal to one. This correspondence can be elaborated upon 
further. In the standard homogeneous isotropic Friedmann-Robertson- Walker model of classical relativistic 
cosmology one uses the line element 

ds 2 = dt 2 - R 2 (t) j-^L + r 2 d9 2 + r 2 sin 2 0# 2 } (5.4) 

with k = 0, ±1 and -ff (to) = {R/ ' R)t Q denoting today's expansion rate as determined from the field equations. 
It is well known that the presence of a small cosmological constant induces an exponential expansion of the 
scale factor at large times. In this very distant future, dominated by a non-vanishing cosmological constant, 

8 In QED the scale dependence of the vacuum polarization effects is of course quite small, with the fine structure constant 
only changing from a(0) ~ 1/137.036 at atomic distances to about a(mz ) ~ 1/128.978 at energies comparable to the Z° mass. 
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an equivalent description can be given in terms of the static isotropic de Sitter metric 

ds 2 = {l - H^r 2 ) dt 2 - (1 - H^r 2 )- 1 dr 2 - r 2 (d9 2 + sin 2 9dcb 2 ) (5.5) 

with a horizon radius H^ = lim^oo H(t). From Einstein's classical field equations one has 

H* 00 =*?£\=± (5.6) 
oo 3 3 \ / 

so the existence of an Hoo is equivalent to assuming the presence of a non- vanishing cosmological constant 
A (here we follow common convention in defining the scaled cosmological constant A, which should not be 
confused with the ultraviolet cutoff). It is presumably this quantity which should be identified with £. Given 
the rather crude nature of our arguments, in the following we shall not distinguish between Ho and H^, 
and simply take Hq 1 ~ 10 28 cm as today's estimate for the size of the visible universe. 9 

The appearance of the renormalization group invariant quantity £ in the quantum evolution of the 
coupling G, a very large quantity by the identification of Eq. (5.3), suggests that the leading scale-dependent 
correction, which gradually increases the strength of the effective gravitational interaction as one goes to 
larger and larger length scales, should be extremely small. One would therefore expect the deviations from 
classical general relativistic behavior for most physical quantities to be in the end practically negligible, at 
least until one reaches very large distances r <~ £. 

At this stage we should comment on an apparent paradox associated with the identification of the 
correlation length £ with 1/H . Naively one would expect, simply on the basis of dimensional arguments, 
that the curvature scale close to the fixed point be determined by the correlation length 10 

K ~ 1/e , (5.7) 

but one cannot in general exclude the appearance of some non-trivial exponent. Indeed one finds for the 
vacuum expectation value of the Ricci scalar (see Eqs. (1.5) and (4.2)) 

W ~ ~~o — 7TTT ( 5 - 8 ) 

with v = 1/3 in four dimensions, and therefore 1Z <~ l/lp £. Only close to two dimensions one recovers the 
classical result, for which v <~ l/(d — 2). 

At first one might be tempted to identify the expectation value of the local scalar curvature with the 
quantity Hq, but further thought reveals that this correspondence is inconsistent with the identification 

9 While the observational evidence for a non-vanishing cosmological constant is quite recent, simplicial lattice theories of 
Euclidean quantum gravity can, as far as one knows, only be formulated with a non-vanishing positive bare cosmological 
constant (A > 0). In the absence of such a constant the path integral does not converge for large edge lengths, and no stable 
ground state exists in the Euclidean theory [8, 9, 16, 13]. 

10 There is at least in principle an even more naive expectation, namely 1Z ~ l/lp, which is excluded though by all numerical 
studies of simplicial lattice gravity. 
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£ = 1/Hq proposed before, and would only be legitimate if the local curvature average 1Z were to indeed 
correctly describe the rotation of vectors, parallel transported around very large loops. But the analogy with 
the local action density < F 2 > in non-Abelian gauge theories seems to suggest that such an identification 
might in fact be incorrect [66], and that the long-distance contribution to the curvature is not given by the 
local average in Eq. (5.8), but should instead be measured directly by computing the parallel transport of 
vectors around very large loops (with characteristic size much larger that the Planck length, A ~ £ 2 ^> 1%), 
but this is laborious and has not been done yet. Indeed at the other end one expects, for very short distances 
comparable to the size of the ultraviolet cutoff, significant fluctuations in the curvature with fluctuations 
of the order of 1Z ~ 1/1%. In other words, the above arguments would suggest that the observable average 
curvature should be scale dependent. 

On the other hand for the curvature correlation at fixed geodesic distance d one obtains from simple 
scaling, and for "short distances" (r <C £), 

<^gR(x)^gR(y)5{\x-y\-d)> c d2( ]_ l/v) ~ £ , (5.9) 

for v = 1/3 in four dimensions, and with A a calculable constant of order one [10]. This last result follows 
almost immediately from the relationship between the curvature-curvature correlation function and the 
second derivative of the partition function with respect to G, which determines the curvature fluctuation 
and thus the curvature correlation function at zero momentum. If one then considers the (scalar) curvature 
R averaged over a very small spherical volume V r = 47rr 3 /3, 

= 7T / d 3 x y/g(x,t) i?(x,t) (5.10) 

one can compute the corresponding variance as 

[S(VgR)} 2 = ^2 f J v d3 y < Vg Vg R(y) >c= |^ • (5.n) 

As a result the r.m.s. fluctuation of y/gR averaged over a small spherical region of size r is given by 

KsfgR) = . (5-12) 

with a Fourier transform power spectrum of the form 

4-TT 2 A 1 

p k = iv5i? k i 2 = -^n • ( 5 - 13 ) 

These results only hold for relatively short distances, and presumably get modified at distances r <~ £. The 
semi-classical answer would look quite different; in this limit v ~ l/(d — 2) and therefore for the curvature 
correlation the distance dependence would tend to 1/d 4 close to d = 2, where it makes sense to make a 
comparison. 
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One can go one step further and use Einstein's field equations to relate the local curvature to the local 
mass density. From the field equations 

R^-lg^R = 8iGT^ (5.14) 

for a perfect fluid 

Tf+u = V9ixv + + P)u^u u (5.15) 
one obtains for the Ricci scalar, in the limit of negligible pressure, 

R{x) w 8ttG p{x) . (5.16) 

As a result one obtains from Eq. (5.9) for the density fluctuations a power law decay of the form 

<p (a 0p( y )> C|x _~ <<£ ^ • (5.17) 

One can list a few other classical general relativistic results which are presumably affected by a running 
gravitational constant. It should be clear from the above discussion that in order for the quantum corrections 
to become quantitatively significant, one needs to look at rather large distance scales comparable to £, or in 
other words r <~ 1/H . In standard classical cosmology one writes 

Ho = ^ [ PA + Pdm + Pb ] (5.18) 

with G usually assumed to be constant, and with the p's representing various density contributions. On the 
l.h.s one usually neglects terms of order k/R^ arising from the curvature of the hypersurface of homogeneity. 
In view of the what has been said before though it seems natural that G(r) in the above expression should 
be taken at the largest length scale Hq 1 . Then one obtains for the overall coefficient a quantity slightly 
larger than the laboratory value ~ G(0), namely G^H^ 1 ) w G(0) (1 + c) > G(0). On the lattice one finds 
a rather small value for c w 0.06. One should recall however, as stated earlier, that while the exponent 
v is universal, the quantity c is not, and in general depends on the specific regularization scheme. More 
specifically, in ordinary lattice gauge theories one finds large but calculable finite renormalization factors, 
relating the lattice gauge coupling to the continuum coupling [67]. A more reasonable expectation would 
therefore be that G(if _1 ) is related to G(0) by a constant of proportionality which is roughly of order one. 
Additional cosmological and astrophysical arguments and proposed tests can be found in [68] . 
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6 Concluding Remarks 



In this paper we have examined some aspects of the connection between lattice and continuum models 
for quantum gravity. In particular the aim of the paper was to elucidate the relationship between the more 
recent simplicial lattice results, which do not rely on the weak field expansion and are therefore inherently 
non-perturbative, and the semiclassical Euclidean functional integral expansion in the continuum. 

The first issue addressed was the very definition of the notion of spin content in the lattice theory 
Proceeding from the Euclidean Feynman path integral approach, we have constructed the lattice analog 
of the semiclassical expansion for the ground state functional of linearized gravity. Two procedures were 
followed. The first procedure relied on constructing directly a lattice expression for the exponent of the 
ground state functional, obtained by transcribing the continuum expression in terms of lattice variables. 
There one proceeds from the lattice expression for the gravitational action, specified on a fixed time slice, 
and supplemented by the appropriate vacuum gauge conditions. A crucial ingredient in this method is 
the correct identification of the correspondence between continuum degrees of freedom (the metric) and 
the lattice variables (the squared edge lengths). The resulting lattice expression is then equivalent to the 
continuum one by construction. 

The second procedure relies instead only on the expression for the lattice gravitational action, as expanded 
in the weak field limit, and determines the explicit lattice form for the ground state functional for linearized 
gravity by performing explicitly the necessary lattice Gaussian functional integrals. The resulting discrete 
expression can then be compared to the continuum one by systematically re-expressing the edge lengths in 
terms of the metric. It is encouraging that the resulting lattice expression completely agrees with what is 
found by using the previous method. 

It is advantageous in performing the above calculation to introduce spin projection operators, which 
separate out the spin zero, spin one and spin two components of the gravitational action. As a by-product 
one can show that the lattice gravitational action only propagates massless spin two (or transverse-traceless) 
degrees of freedom in the weak field limit, as is the case in the continuum. Furthermore, as expected the 
lattice ground state functional for linearized gravity only contains these physical modes. 

The explicit construction of the ground state wave functional of linearized lattice gravity in terms of 
the lattice transverse-traceless modes makes it possible at least in principle to compare the lattice and 
continuum results in the limit of small curvatures. After imposing appropriate boundary conditions at 
infinity by suitably restricting the values for the edge lengths on the lattice boundary such that the deficit 
angle is zero there, one would then have to enforce as well the lattice vacuum gauge conditions of Eq. (2.48) so 
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as to make contact with the semiclassical lattice functional of Eq. (2.42). Since no gauge fixing is required for 
determining invariant averages obtained via the partition function of Eq. (1.1), the gauge conditions would 
have to be imposed on each edge length configuration, by progressively applying local gauge transformations 
[14]. But it is expected that after such a transformation the edge distributions on a fixed time slice should 
follow closely the distribution of Eq. (2.42), if indeed as expected the only surviving physical modes are 
transverse traceless. 

In subsequent sections of the paper we have systematically examined the relationship between recent 
non-perturbative results obtained in the lattice theory and the corresponding calculations as performed in 
the continuum theory. The latter suggest the presence of a non-trivial ultraviolet fixed point in G, and in 
some cases have led to definite predictions for the universal critical exponent of quantum gravitation, which 
can therefore be compared - even quantitatively - to the lattice results. 

Besides relying on the recent lattice and continuum results for quantum gravitation, one can also in- 
dependently try to estimate the gravitational scaling dimensions based on what is known based on exact 
and approximate renormalization group methods for spin zero (self-interacting scalar field in four dimen- 
sions) and spin one ( Abelian non-compact gauge theories) , for which a wealth of information is available on 
the critical indices. We have argued that these results too are remarkably consistent with what is known 
about the gravitational exponents in four dimensions, both from the discrete as well as from the continuum 
side. We have also presented a simple geometric argument which interprets the value for the gravitational 
exponent v~ x = 3. 

In the last section of the paper we have discussed some (almost immediate) physical implications of 
recent lattice and continuum results, with an emphasis on the small expected deviations from classical 
general relativity expected at sufficiently large scales due to the running of G. We have argued that it is an 
almost inevitable consequence of the existence of an ultraviolet fixed point that the gravitational coupling 
becomes scale dependent, with power law corrections involving the correlation length. In analogy with 
non- Abelian gauge theories, and in the absence of any other likely physical candidate, it seem natural to 
identify the non-perturbative scale determining the size of deviations from classical gravitation with 1/H , 
as suggested in [15]. 
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